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Abstract 

Multiphoton femtosecond coherent control of is used for implementing innovative photo-induced 
analog coherent computation that generally might be a basis for future "smart hardware". The 
specific implemented computational task the classification of an unknown sequence into one of the 
three groups: (i) a constant sequence that is composed of identical numbers, (ii) a sequence that is 
antisymmetric around a given point, or (iii) neither. The input sequence is encoded into the spec- 
tral phases of a broadband femtosecond pulse and the computational task is being carried out by 
the multiphoton nonlinear response of the irradiated physical system. Here, it is the simultaneous 
coherent two- and three-photon absorption in atomic sodium (Na). The corresponding computa- 
tional resources are the manifold of initial-to-final multiphoton excitation pathways photo-induced 
by the broad spectrum of the irradiating femtosecond pulse. The answer is obtained by measur- 
ing only two observables, which are two state-populations excited via the two- and three-photon 
absorption processes. Hence, the computational task is accomplished in a constant number of op- 
erations irrespective of the sequence length. As such, the presented scheme, where the femtosecond 
pulse serves as a query and the irradiated physical system serves as an oracle, provides a sufficient 
gain in the computational complexity. 

PACS numbers: 32.80.Qk, 42.50.Ex 
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Analog quantum computation (AQC) utilizes the similarity between the formal math- 
ematical description of the physical phenomena and the formulation of the problem being 
solved. The information is processed by exploiting the physical phenomena directly: it is the 
system evolution that processes the information and performs a computational task 
AQC is aimed to solve otherwise intractable problems, once the similarity is found between 
the mutual relations of the problem parameters and the corresponding physical quantities 
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gL^, y, 9, Q, lll(. Here we exploit the femtosecond coherent con- 



17( 1 of multiphoton absorption to perform a specific computational 



task: classification of the unknown sequence in to the three groups: a sequence composed 
of the identical components, a sequence that is antisymmetric about a specific point, or 
neither. The concept of using the femtosecond pulse as a query and an excited system as 
an oracle provides a sufficient gain in the computational complexity. The experimentally 
demonstrated approach requires only two measurements to accomplish the computational 
task irrespectively of the sequence length. 

Following from their coherence over a broadband spectrum, femtosecond pulses are used 



as a unique tool for coherent control and coherent spectroscopy of matter 16|, 
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22|. The broadband spectrum provides a variety of photoinduced pathways that are 



used for the multiphoton excitation. All of the such pathways interfere between them and 
contribute to the finally excited amplitude: constructive interference enhances the excited 



amplitude, whi 
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e destructive interference suppresses it. The recently developed pulse-shaping 
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25l ] enables access to the parameters of different frequency components of 



the pulse independently. Among the parameters that could be adjusted are phase, amplitude, 
and polarization of the frequency components. Control over these parameters is directly 
translated to the control over the interference between the photoinduced pathways and, 
hence, to the control over the excited amplitude. 

The basic idea of the presented here experiments is to use the interference pathways 
of the multiphoton excitation as the computational and information processing tool. The 
unknown information to be processed is translated into the transformation f/=exp[i$(co')] 
which is applied to the electric field of the excitation femtosecond pulse with spectral am- 
plitude This transformation results in a shaped pulse |£'(to')|exp[i$(c<j)], where ^{u) 
is a spectral phase of frequency uj. The phase associated with each multiphoton pathway is 
equal to the sum of the phases of all the involved photons. All the possible combinations 
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of the excitation pathways provided by the shaped pulse interfere together during the exci- 
tation, i.e. evaluated at a single event, as opposed to the classical information processing, 
where each information component is evaluated separately. After the excitation, the re- 
sulted excited amplitude contains the information on the nature of the interference. Relying 
on the knowledge of the light-matter interaction mechanism, measurement of the excited 
amplitude can characterize the information originally used to shape the excitation pulse. 
Different physical systems respond in a different way to the excitation with a given pulse 
shape. This is why the system used for the processing should be carefully chosen to suit 
the specific computational task. It is a general and very important attribute of AQC - its 
inherent problem-specific nature: one should always find a physical system with quantities 
which mutual relation corresponds to these of the problem to solve. 

The computational task that we consider is a classification of the unknown sequence in 
to the three groups: a sequence composed of the identical components, which is referred 
as constant sequence, a sequence that is antisymmetric about a specific point - antisym- 
metric sequence, or neither, referred as random. For that purpose we encode the subjected 
sequences into the phase patterns ^const, ^antisym, and ^random of the broadband spectrum. 
Each component of the sequence is translated into the phase magnitude of the corresponding 
frequency using the pulse shaping technique. The physical process that suits the compu- 
tational needs is the two-channel femtosecond multiphoton absorption. In the proposed 
scenario the amplitudes of the excitation channels characterize the type of the sequence 
encoded into the transformation U. The different response by the different channels to the 
applied transformation is a key property of the system that solves the presented classification 
task. 

One of the excitation channels exploited for our computational task is a non-resonant two- 
photon absorption channel Q- The excited amplitude in this case is proportional to 
A^'^\n) = E{^l/2 — a)E(^l/2 + a)da. It interferes all the possible two-photon pathways 
with transition frequency fl, each such pathway is composed of two photons with frequencies 
Q/2 — a and Q/2 + a, i.e., two frequencies located symmetrically around Q/2. The phase 
associated with each pathway is $(^2/2 — a) + $(i7/2 + a). For the transform-limited (TL) 
pulse, which is the shortest pulse for a given spectrum, $(to') = and all the excitation 
pathways interfere constructively. The same applies to the transformations corresponding 
to the constant sequences f/™"^* with $(ti;)=const as they add only a global phase factor. 
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identical to all the interfering pathways. The two-photon non-resonant channel is also in- 
variant to the transformations corresponding to the antisymmetric sequences f/""-*«'^y™- ^ while 
the point of the inversion corresponds to half the two-photon transition frequency VL/2. For 
that case $(^2/2 — a) = — $(i7/2 + a) holds, and all the available pathways obtain zero 
relative phase. Hence, the application of Jj^-^^^^v"^ implies constructive interference between 
all the pathways. On the other hand, along any other possible channel without special sen- 
sitivity to the symmetry of the transformation, the induced excitation pathways interfere 
destructively. Such channel is not invariant to (/"■"■^^■'^y^^ but to f/'^°"''*. 

Combining all the described here together, we build a truth table for the sequence classi- 
fication: constant sequences induce fully constructive interference along both channels and 
they are associated with a maximal amplitude of both excitation channels, which is actu- 
ally the same as for the TL pulse. The antisymmetric sequences are associated with fully 
constructive interference along the two-photon non-resonant channel, and with destructive 
interference along the channel without any symmetrical properties. The random sequences 
induce destructive interference along both channels and the excited amplitude by the corre- 
sponding shaped pulses is sufficiently less then by the TL pulse. To demonstrate this idea 
experimentally, we have chosen atomic two-channel system, which includes non-resonant 
two-photon absorption channel coherently incorporated within a resonant-mediated (2+1) 
three photon absorption channel. 

The two-channel excitation considered here is shown schematically in Fig. [TJ It involves 
an initial ground state \g) and two excited states and I/2). The \g) and states are 
coupled by a non-resonant two-photon coupling provided by a manifold of states \n) that 
are far from resonance, while the and I/2) states are coupled resonantly by one-photon 
coupling. Hence, when irradiated with a weak (shaped) femtosecond pulse of temporal 
electric field e{t), two multiphoton absorption channels are induced simultaneously: a non- 
resonant two-photon absorption channel from \g) to is induced simultaneously with a 
resonance- mediated (2+1) three-photon absorption channel from l^f) to I/2) with as 
an intermediate state. The amplitude of the latter channel is determined by both intra- 
group and inter-group interferences corresponding to on-resonant and near-resonant with 
1/1) state three-photon excitation pathways respectively. The pathways induced by the pulse 
shaped with jj^-'^^'^^v'^ interfere destructively along this channel. Although the TL pulse and, 
hence, the shaped pulses corresponding to f/™"^* do not induce the maximal constructive 
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interference for the (2+1) resonance-mediated channel, the pulse should be shaped in a very 
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special way to enhance the amplitude excited by the TL pulse |20|, 
case of the present experiment. 

Alternatively, another two-photon non-resonant channel can be chosen for that purpose 
with different transition frequency ^Ichi 7^ ^ch2- Obviously, while the channel 1 is sensitive 
to the sequences antisymmetric about Qchif^, the channel 2 is sensitive to the antisymmetric 
sequence with the inversion point at Qch-^f^- The corresponding physical system is a second 
harmonic crystal, with each frequency of the produced harmonics being a different excitation 
channel measured separately. The choice of such a crystal provides a wide tunability upon 
the point of antisymmetric inversion of the sequences to be checked, while the sensitivity 
to this specific point (the precision of the classification) in this case is set by the resolution 
of the detection system. In the atomic case, this sensitivity is extremely high due to the 
narrow width of the atomic lines. 

Suppose now that unknown sequence f{x) is optically encoded into the phase pattern 
$(co') of the pulse using a pulse shaping technique. We check whether there is an inversion 
point Xq corresponding to the frequency of the non-resonant two-photon transition Q/2 = 
The processing of the information is done upon the excitation of the atomic system 
with the shaped pulse from its ground state \g) to the excited states and I/2). The 
population of the excited states Pj^ and Pj^ is used as a readout channel and compared with 
the population excited by the TL pulse. If Pf^j{x) = Pfi,TL and Pf2,f{x) = Pf2,TL, the pulse 
spectral phase is a constant pattern which represents encoded sequence /(x) = const. If 
PfiJix) ~ Pfi^TL, but Pf^ji^x) 7^ Pf2,TL, the pulse spectral phase is antisymmetric about Ufg 
the encoded sequence satisfies /(xq — a) = — /(xq + a). If Pfij{x) 7^ Pfi,TL and Pf^j(^x) 7^ 
Pf2,TLi the corresponding phase pattern is neither antisymmetric nor constant and is referred 
here as random. Overall, the computational task is performed by comparing the outcome of 
the two specific measurements only, irrespectively of the length of sequence to be classified. 

The physical model system of the study is the sodium (Na) atom [261], with the 3s 
ground state as \g), the 4s state as and the 7p state as I/2) (see Fig. [1]). The transition 
frequency ijJf^,g = ^jJiis,zs = 25740 cm^^ corresponds to two 777- nm photons and the transition 
frequency t^/a./i = ^7p,is = 12801 cm~^ corresponds to a 781.2-nm photon. The 3s-4s non- 
resonant two-photon coupling originates from the manifold of p-states, particularly from the 
3p state [c<j3p 3s~16978 cm^^ (589 nm)]. The sodium is irradiated with phase-shaped linearly- 



5 



polarized femtosecond pulses having a Gaussian intensity spectrum centered around 780 nm 
(12821 cm~^) with 5.8-nm (95-cm~^) bandwidth (~180-fs TL duration). 

Experimentally, a sodium vapor in a heated cell is irradiated with such laser pulses, after 
they undergo shaping in an optical setup incorporating a pixelated liquid-crystal spatial 
light phase modulator [23|. The effective spectral shaping resolution is 6uJshaping='^-05 cm~^ 
(0.125 nm) per pixel. The peak intensity of the TL pulse is below 10^ W/cm^. Following 
the interaction with a pulse, the population excited to the 4s state radiatively decays to the 
lower 3p state, which then decays to the 3s ground state. The 3p-3s fluorescence serves as 
the relative measure for the total 4s population P/^ = P^s- The population excited to the 
7p state undergoes radiative and collisional decay to lower excited states, including the M, 
5d, 6d, and 6s states. The fluorescence emitted in their decay to the 3p state serves as 
the relative measure for the total 7p population Pf^ = Pip. The fluorescence is measured 
using a spectrometer coupled to a time-gated camera system. The 3p-3s fluorescence part 
originating from the 4s state is discriminated from the part originating from the 7p state by 
using a proper detection gate width, utilizing the different time scales of the 4s-to-3p and 
7p-to-3p decays. 

Fig. [2] presents experimental results measured for the randomly generated 20,000 se- 
quences encoded into the phase patterns. Half of the sequences are urged to be antisymmet- 
ric about uj4s,3s- Each pattern is composed of 26 bins which can obtain values randomized 
in the [— 7r/2,7r/2] interval. On Fig. [2] panel (a) presented a two-dimensional histogram of 
the measured population: x-axis represents the population of 4s state, y-axis represents the 
population of 7p state. The color code and the topological curves represent the occurrence 
of phase patterns with the corresponding signal. All the results presented here are given 
relative to the absorption induced by the transform-limited (TL) pulse. For convenience, we 
introduce the TL- normalized absorption measures P/^ = P/^/P/^^tl and P/j = P/2/-P/2,rL- 
The three regions, corresponding to the constant, antisymmetric, and random phase patterns 
are clearly distinguished on the histogram. The region with the peak at (P4s=l, P7p=0.5) 
correspond to the antisymmetric sequences, encoded to the phase patterns, while the region 
with the peak at (P4s=0.25, P7p=0.25) correspond to the random sequences. The region 
with the peak at (P4s=l, P7p=l) correspond to the constant sequences with the populations 
identical to the TL pulse. The broadening of that peak directly reflects the signal-to-noise 
experimental ratio. As seen from the figure, the measured population, corresponding to the 
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sequence encoded to the phase pattern of the excitation pulse can be used for the classifica- 
tion of the sequence by its properties, based on the presented above truth table. 

On Fig. [2] panels (b) and (c) we draw the projections of the two-dimensional histogram 
of panel (a) on the 7p and 4s population axis. As seen from panel (c), the 4s population 
signal clearly distinguish between the random sequences and the rest: 100% of the random 
sequences produce signal less then the corresponding TL pulse. However, this signal do 
not distinguish between the constant and antisymmetric sequences. The latter is done by 
analyzing the 7p population signal [see panel (b)]. For the detailed analysis of the distri- 
bution we plot the cumulative sum function corresponding to the antisymmetric sequences 
in dashed black line and obtain that 96% of the sequences produces 7p population signal 
less then the corresponding TL pulse. It means that the experimental false identification 
of antisymmetric sequence as a constant one or vise versa is below 4%. This value can be 
improved by improving the signal-to-noise ratio. 

We have checked the sensitivity of the method by shifting the phase pattern center of 
symmetry by one pixel apart from uj4s,3s- The corresponding results are shown on Fig. [2]panel 
(b). The distribution of the antisymmetric patterns is moved from the region of P^s = 1 
to the region corresponding to random patterns. The reason for the high sensitivity of the 
signal to the exact location of the inversion point is in the very narrow nature of the atomic 
resonance. While the phase pattern is antisymmetric about ujis,3s/'^ it induces completely 
constructive interference of all the photoinduced pathways to the |4s) state. When the 
inversion point is shifted by A, the maximally constructive interference would be induced 
for the state corresponding to the frequency uj4s,3s + A. In our experiment we have shifted 
the pattern by A = 6u!shaping='^ -^^cm^^ . This shift was sufficient to reduce the constructive 
interference resulted in the reduced signal measured for that case. 

To summarize, we have demonstrated an efficient solution to the unknown sequence 
classification problem. The classification is into three possible groups: a constant sequence, 
a sequence that is antisymmetric about specific point, and a random sequence. The sequences 
are encoded into the phases of the broadband femtosecond pulses. In the spirit of the analog 
quantum computation, the photoinduced dynamics of the excited atomic sodium is used 
to perform the computational task. The solution to the problem is obtained within the 
measurement of the two excitation channels irrespectively of the sequence length. 
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FIG. 1: The two-channel femtosecond excitation scheme of Na, including a non-resonant two- 
photon absorption channel from \g)='is to \fi)=As and a resonance-mediated (2+1) three-photon 
absorption channel from \g)='is to |/2)=7p via |/i)=4s. Shown in solid lines are examples of two- 
and three-photon pathways. The latter are either on resonance or near resonance with (with 
detuning 5). The decay to the lower excited states is shown by dotted lines, the measured radiative 
decay indicating the population of the 7p and 4s states is shown by dashed lines. 
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FIG. 2: Experimental results for the 20,000 sequences encoded to the phase patterns of the shaped 
pulses. Each pattern has 26 bins with the values randomly generated within the interval of [— vr, vr]. 
Half of the sequences were urged to be antisymmetric about uj/^s,3s- Panel (a): Presented is two- 
dimensional histogram of the normalized to the unshaped (TL) pulse populations of 4s and 7p 
states. Three regions corresponding to the constant, antisymmetric, and random sequences are 
clearly distinguished. Panel (b): Projections of the panel (a) histogram on the axis of 7p population. 
Dashed black line indicates the cumulative sum of antisymmetric sequences histogram. At (x = 1) 
it obtains the value of 0.96. Panel (c): Projections of the panel (a) histogram on the 4s population 
axis. Panel (d): Two-dimensional histogram of antisymmetric sequences shifted by A = Slo shaping 
from the original center of antisymmetry u;4s_3s/2. The resulted sufficient shift of the distribution 
indicates the high sensitivity of the proposed method. 
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